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\^ , We study the possibility to detect the gravitational wave background generated by all the neutron 

stars in the Galaxy with only one gravitational wave interferometric detector. The proposed strategy 
ON ■ consists in squaring the detector's output and searching for a sidereal modulation. The shape of 

the squared signal is computed for a disk and a halo distribution of neutron stars. The required 
■ noise stability of the interferometric detector is discussed. We argue that a possible population of 

old neutron stars, originating from a high stellar formation rate at the birth of the Galaxy and not 
emitting as radio pulsars, could be detected by the proposed technique in the low frequency range 
of interferometric experiments. 

m 

(N ! PACS number(s): 04.30.Db,97.60Jd,97.60.Gb 



I. INTRODUCTION 



Rotating neutron stars (NS) are possible astrophysical sources of gravitational radiation in the frequency range 
of the interferometric detectors LIGO, VIRGO and GEO600 currently under construction [|J-B[. Indeed, provided it 
deviates from axisymmetry, a rotating NS emits continuous wave (CW) gravitational radiation mainly at its rotation 
frequency and twice this frequency g. The non-axisymmetric shape of a NS can be caused by anisotropic stresses 
• from the nuclear interactions, irregularities in the solid crust ("mountains") [^|J|], the internal magnetic field 
Oh| some precessional motion ||, or the development of triaxial instabilities (for the most rapidly rotating NS): the 
Chandrasekhar-Friedman-Schutz instability Jl^;|ll| and the MacLaurin-Jacobi type instability induced by viscosity 
JHqJIil . Estimates of upper bounds on the individual gravitational wave amplitude from a sample of 334 observed 
NS (radio pulsars) have been provided by Barone et al. JTJ] (see ref. (resp. ref. for recent values based on 
a sample of 558 (resp. 706) pulsars). The amplitude of gravitational waves (GW) emitted by a rotating NS can be 
expressed in terms of the NS rotation period P, its distanc e to the Earth r, its moment of inertia I about the rotation 
axis and its ellipticity (triaxial deformation) e as (cf. Eq. JA5| ) below) 



h = 4.21 x 10- 24 PfWf *)('). 

VP/ V r J \.10 38 kgm 2 y V 10- 6 J 



(1) 



The crucial parameter entering this formula is the ellipticity e. Its value depends on the physical mechanism that 
makes the star non-axisymmetric (cf. the references given above) and is highly uncertain. Upper bounds on e can 
however be derived from the observed slowing down (P) of pulsars by assuming that this latter is entirely due to 
the loss of angular momentum by gravitational radiation. Let us note that this provides an absolute upper bound; 
most of the P is usually thought to result instead from losses via electromagnetic radiation and/or magnetosphcric 
acceleration of charged particles — at least for Crab-like pulsars. The maximum values of e obtained in this way, 
as well as the corresponding maximum values of ho, are given in Table § for five pulsars. The first three of them 
correspond to the three highest values of /lo.max among the 706 pulsars of the catalog by Taylor et al. |l7],[l8|]. The 
two remaining entries are two millisecond pulsars: the second fastest one, PSR B1957+20, and the nearby pulsar PSR 
J0437-4715. The figures ho, max ~ 10 -24 for the Crab-like pulsars (three first entries in Table ||) are almost certainly 
too optimistic because, as already said, electromagnetic phenomena can be invoked to explain most of, if not all, the 
observed pulsar spin-down. Besides, one may notice, following New et al. ITJ], that if the mean ellipticity of pulsars 
is taken to be of the order of the e max of millisecond pulsars, i.e. e ~ 10~ 9 (cf. Table |), then the Crab pulsar reveals 
to be a much worse candidate than PSR B1957+20, as it can be seen by setting e = 1CP 9 in Eq. ([!]) for these objects: 

h c r * b _ 2 x 1Q -30 vergus h ^7+20 _ x x 10 -27_ 
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The detectability of individual NS by existin g an d future gravitational wave detectors has been discussed b y v arious 
authors, including Schutz f2C(] ; Jotania et al. |j2l| , ^2| . Suzuki 0], New et al. E3], and Dhurandhar et al. |24| . For 
VIRGO-like instruments, it can be hoped that any pulsar that produces a GW amplitude on Earth, ho, greater than 
10 -26 in the frequency bandwidth where the sensitivity of the detector is better than 1CP 22 Hz -1 / 2 , can be detected 
with three years of integration |l6[ . The last column of Table [j] gives the minimum value of e required to produce 
ho > 10~ 26 . Notice that for Crab-like pulsars this value is below 1% of the maximum ellipticity allowed by the 
measured spin-down rate of the pulsar, whereas for millisecond pulsars both values are of the same order. 

In the present article we study the possibility to detect the total CW emission from the whole population of NS in 
our Galaxy by using only one LIGO/ VIRGO type detector. The proposed strategy consists in measuring the square 
of the gravitational signal, h 2 , and in detecting the sidereal modulation of this signal which results from the directivity 
of the detector and the anisotropy of the NS distribution. This technique (referred hereafter as quadratic detection) is 
very similar to the one used in radioastronomy when only one antenna is used and differs from the strategy proposed 
by Schutz EfJ that consists in searching for NS one by one within a 4-dimensional space (frequency, phase and position 
on the sky) (technique referred hereafter as linear detection) . The advantages and drawbacks of the two strategies are 
discussed and it will be shown that if the number of the emitting stars is larger than 2 x 10 6 , then the signal squaring 
technique is more convenient and more computer time saving than the linear one. Actually, the two techniques appear 
to be complementary. 

The paper is organized as follows: in Sect. [0] the statistical properties of the squared signal are briefly recalled. The 
efficiency of the method as a function of the NS number and of the range of frequency at which they are supposed 
to radiate is computed, the comparison with the linear analysis for single NS search being performed in Sect. III. 
In Sect. [V, the constraints on the stability of the instrumental noise are computed. In Sect. we try to evaluate 
the number of gravitationally emitting NS in our Galaxy and the range of frequencies at which they are supposed to 
radiate. Sect. VI contains the conclusions. 



II. STATISTICAL PROPERTIES OF THE SQUARED SIGNAL 

The response h(t) of an interferometric detector to the gravitational radiation emitted by N NS spread out on the 
sky is detailed in Appendix |A|. For the purpose of this section, let us consider that hit) is the sum of N ele ment ary 
periodic functions hi with unknown frequencies Vi and phases <fii, the precise form being given by Eqs. ( |A2| )-( [A5| ). 

The time average value of h(t), (h(t)), is zero, but the average of its square is not. Therefore the strategy to 
detect the gravitational emission of an ensemble of NS, the frequencies of which spread out from v\ to v-i, consists in 
measuring the square of the s igna l h 2 (t). It is easy to see that (h 2 (t)) is proportional to the sum of the squares of the 
shear from each NS (cf. Eq. flA13p): 



(h 2 (t)) 



N 

£ 

i=l 



Cti(t)'- 



(2) 



where ri is the distance of the i-th NS, Ai its gravitational wave amplitude at one unit distance, OLi(t) some factor 
involving the direction of the NS and the polarization of its radiation with respect to the detector arms. For the 
purpose of the present discussion, Eq. (|J) can be recast in the following approximate form 



(h 2 (t))^K(t)N — 



(3) 



where A 2 is the mean value of A 2 , 
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is the inverse-square average distance of the NS population and K(t) is a time varying factor, with the periodicity 
of one sidereal day. Its non-constancy is induced by the directivity of the detector and the anisotropy of the NS 
distribution. 

The study of the efficiency of the quadratic technique is made easy by the close analogy with the radioastronomy 
observation technique. Radioastronomers measure the square of the electromagnetic field emitted by a large ensemble 
of collective modes of a plasma. They are confronted to the problem of extracting some signal from the noise of the 
receiver. One possible technique consists in scanning the sky around the source, and in measuring the differences of 
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the total noise on/off of the source. In the case of gravitational waves, the source is scanned by the interferometric 
detector via the rotation of the Earth. 

Let us also notice that the problem of detection of the gravitational wave emission from an ensemble of NS is very 
close to the search for a gravitational cosmological background discussed in great details by Flanagan (see also 
p6|). The main and basic difference is the daily modulation of the signal from the NS ensemble, which allows the 
search to be performed with only one detector instead of two detectors required for the cosmological background 
detection p5|] . 

Let s(t) = h(t) + n(t) be the output of the detector, hit) being the gravitational radiation signal and n(t) the 
detector's noise. The frequency bandwidth of the detector is supposed to be Av = v<i — v\. By squaring and by 
averaging s(t) we obtain: 

(s 2 (t)) = {h 2 (t))+{n 2 (t)) (5) 

because noise and signal are not correlated and consequently the average of the cross product (2n(t)h(t)) vanishes. 
Here the average must be understood as the average on the outputs of an infinite number of detectors. 
If the noise is stationary, (n 2 (t)) is a constant, and its value is 

(n 2 (t)) = / G(y)h 2 (y)dv (6) 

J V\ 

where h 2 (v) and G(v) are respectively the noise power per Hz and the frequency response of the detector. For 
simplicity we shall assume the noise being white (n independent of v and G(y) = 1 for v\ < v < v<i ). The above 
expression then reads 

(n 2 (t))=h 2 Av (7) 

The quantity (h 2 (t)) is time dependent: there are high frequency time variations with a typical frequency of NS 
rotation, and a low frequency time variation due to the slow change of the orientation of the interferometer arms with 
respect to the Galaxy induced by the Earth rotation. This latter frequency, ^ S id, is equal to the inverse of one sidereal 
day: = 86164.092055 s. It is this modulation that must be searched for. 

In practice, observations are performed with only one detector; consequently the ensemble average () must be 
replaced by an average on time. In what follows we shall consider a simplified case in order to allow the reader who 
is not familiar with the radioastronomy technique, to understand the basic ideas. The reader will find more details 
in the excellent book by Kraus |^7j and in the already quoted paper by Flanagan ^5|. Let T be the averaging time. 
We can write 



t+T n 2 {t')dt' = {n 2 {t))+^(t) , (8) 



where tp(t) is a random function. Under the ergodicity hypothesis ip(t) vanishes when T — > oo. When T is finite, ip{t) 
has the following statistical properties: 

M*)>=0, <^ 2 (t)) = Cn i Av/T (9) 

Here C is a constant of the order of unity and depends on G(v). For the rectangular filter considered above C = 2. 
Values of C for different filters can be found in Chapter 7 of ref . p7] . 

If T is much shorter than one sidereal day, (ft, 2 (i)) can be considered as constant, and the signal-to-noise ratio i? qua d 
can be easily computed from Eqs. (||) and (||): 

NA 2 NA 2 VT , s 

qUad ~ D\JWWl ~ D l n 2 VCA^ 

If H = N/Af is the number of the NS per unit frequency, i? qua d reads 



HA 2 VAvT , , 

Rquad ~ n 2 D 2 VC 



Note that i? qua d is magnified by the factor v Ai/T well known by radioastronomers. 
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If T is much longer than one sidereal day, we have to take into account the periodicity of the signal. Let us suppose 
for simplicity^] that the factor K(t) appearing in Eq. (0) is harmonic with period P s id equal to one sidereal day and 
with a known phase: (h 2 (t)) = N A 2 cos(2irt / P^d) / D 2 . In this case the best way to proceed is to make a Fourier 
transform of the squared signal. The signal-to-noise ratio i? qua d is then: 

NA 2 

-Rquad — , (12) 

where 5v = 2/T is the width of a single bin of the Fourier transform and il> 2 (v) is the spectral power of the noise 
ijj(t) 2 . At very low frequency V> 2 (y) = Cn 4 Ai/, so that Eq. (O) becomes 



^ uad = N ^S ^ ' (13) 

In the realistic case (cf. AppendixJAj) , the signal is periodic (period P s id) but not harmonic. Its shape depends on the 
Galactic NS distribution. Figures [ljand § show its variation during one sidereal day for two different NS distributions. 
These shapes can be used as templates to optimize the extraction of the signal. 

III. COMPARISON WITH THE LINEAR SEARCH FOR A SINGLE NS 

Let us compare the efficiency of the quadratic analysis with respect to the linear one proposed by Schutz |20| for 
single NS searches. If the frequency and the position of the NS one searches for is known, the signal-to-noise ratio of 
the linear technique is given by 



where D; is the distance to the individual NS that is searched for. The factor 4 in the denominator is due to the 
product of the bandwidth times an extra-factor 2 coming from the fact that the phase is not known. 

The ratio E = R\i n / W -R qU ad is a good quantity to characterize the advantages and drawbacks of the two methods. 
We have, from Eqs. (jl|) and @, 

D q (CAuT)^ 

~D[ 23/4 JVV2 ■ ( 15 ) 

Taking Av = 1 kHz, T = 1 yr, C = 2 and Dj = D q , the quadratic technique appears to be more efficient (E < 1) than 
the linear one if the number of NS is larger than 9 x 10 4 . The above result is based on the two underlying hypothesis: 
(i) the frequency and the position of the isolated NS are known exactly, (ii) the distance D t of the isolated NS is equal 
to the averaged distance D q = {l/r 2 )^ 1 / 2 of the NS of the ensemble. 

Let us consider now the case D q ^ D{. D„ depends on the distribution of NS in the Galaxy. If their distribution 
corresponds to the Galactic disk (see Sect. |v|), then D q — 5.1 kpc (value resulting from the distribution function 
( A21| ) with the parameters Rq = 3.8 kpc and zq = 0.5 kpc). D\ can be taken to be the distance of the closest NS. For 



example, for the nearby millisecond pulsar, PSR J0437-4715 J28|, D\ = 140 pc. Let us take D\ = 100 pc. The two 
methods become then equivalent for N = 2.2 x 10 s NS. 

However, in searching for a single NS, its position has to be known with an accuracy high enough to compensate for 
the Doppler shift of the frequency induced by the rotation of the Earth and its motion around the Sun. The last effect 
is the most important one. The accuracy of the declination 5 and in the azimuth <f> is about 6<f> = 0i?0t' max c where 
R® ! "max and c are respectively the radius of the Earth orbit (150 x 10 6 km), the expected maximum frequency of 
the NS and the speed of light. A similar precision is needed on the declination of the source. This means that the 
sky must divided in about 4 x 10 10 (v/100 Hz) 2 boxes and we have to try to detect a periodic source (by Fourier 
transform) in each box, after compensation of the variation of the frequency of the source due to the Earth motion. 
In the above rough analysis, we have neglected the less important Doppler shift induced by the Earth rotation. In 



1 The actual time variation of (h 2 (t)) is given in Appendix |a|. 
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this section, we do not discuss the technical possibility of performing more than 10 10 Fourier transforms each of them 
containing 3 x 10 9 (i/ max /100 Hz)(T/l yr) bins. The detection will be considered as positive if the probability of a 
random fluctuation of the signal in all the boxes times the number of bins is less than a prefixed value. If we take the 
probability to be lower than .16 (corresponding to the 1 a criterion), the value of the corresponding signal-to-noise 
ratio -Run is about 9.5 for v max — 100 Hz and 10.5 for f max = 1 kHz. Consequently Eq. jl5| ) becomes 

The two methods become then equivalent for TV = 2.2 x 10 6 . 



IV. STABILITY OF THE NOISE 



A. Non-stationary noise 



The above results hold under the hypothesis that the term (n 2 (t)) in Eq. (Q) is constant, i.e. that the noise is 
stationary. Now the actual noise is not stationary: low frequency fluctuations are always present. Let us recall that 
the sources of noise in interfcrometric detectors are the photon shot noise (at high frequency) and the Brownian noise 
of the mirrors (at low frequency). There are at least two sources of low frequency fluctuations: (i) the fluctuations of 
the optical power of the laser and (ii) the temperature fluctuations of the mirrors. When the interferometer is in lock, 
the fluctuations of the laser power can be taken under control within few 10 -5 and are therefore not dangerous at 
all. However, the interferometer may be falling out of lock occasionally or the laser may be shut down and switched 
on some time later. In either case, this will induce a temperature fluctuation in the mirrors. However its amplitude 
will be at most 0.1 K p9|. This falls within the range of the required accuracy on the temperature measurement of 



non-periodic fluctuations, as derived in Sect. IV B below. 

We shall thus consider only the mirror temperature fluctuations and estimate the constraints they impose. We will 
discuss the possibility of monitoring the temperature fluctuation to keep their effects under control. 

In what follows, we suppose that the fluctuations in the term (n 2 (t)) of Eq. @ have a small amplitude and that 
their typical time scale is much longer than Let us introduced the new random variable a(t) by 



(n 2 (t)) = (n 2 (t))(l + a(t)) (17) 

where the symbol X means the average on time (if it exists). The hypothesis of small amplitude fluctuations implies 
that yj (a 2 (t)) < 1. As stated above, we consider only the contribution to the noise due to the temperature. In this 
case, the noise n 2 {v) is proportional to the temperature 6 (|(]] ||l[] : 

h 2 {v)=f3Q. (18) 

Under the above hypotheses, Eq. (||) reads: 

1 £ +T n 2 (t>) dt' = WW) (1 + «(*)) + m (19) 
where tp(t) has the same properties as that stated in Eq. (||). Within quadratic terms in a (t), Eq. @ reads 

i r t+T 

- / n 2 {t') dt 1 = h 2 Ai> + n 2 Av a(t) + ip(t) . (20) 

-L .1 + 



The random variable a(t) has a zero mean value, a(t) — 0, and its variance a 2 (t) is proportional to the temperature 
fluctuations of the mirrors. Consequently the typical time scale for a(t) is of the order of one hour, much longer 
than 1/fi ~ 0.1 s. The typical amplitude of the temperature variation is a few kelvins, much smaller than the room 
temperature (300 K), therefore the above hypotheses are fulfilled. It is worth to note that in defining a(t) and its 
statistical properties we have averaged on time, instead of averaging on an infinite ensemble of identical detectors. The 
reason is that for identical detectors in the same environment, the daily variation of the temperature, and consequently 
the amplitude of the thermal noise, is the same for all the detectors. 
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B. Non-periodic temperature fluctuations 



The fluctuations of the mirrors temperature are not known today. Only in situ measurements of the temperature, 
once the detector will be operational, will provide us with the statistical properties of a(t). However, let us try to 
guess the temperature trends of the mirrors, in order to see if the precision needed on the temperature control can 
be achieved with the present technology. One reasonable hypothesis is to take a power law for the Fourier spectrum 
of the temperature fluctuations <50 (v) 2 : 

2 / 

a I v 



SQ(uY = — — (21) 

where a, 7 and vt are some constants. These parameter can be estimated in the following way. According to Eq. (|2 
the temperature variation on a time scale r is given by 

/ \ — 7 



((5e 2 > = — / — dv> = — , (22) 

*T J v \vi J 7 - 1 \vt J 

where v = 1/r. As already said, we do not have yet any real measurements of the temperature variations of the 
detector's mirrors. Consequently we shall proceed by guessing the temperature variations in the environment of the 
experiment. It seems reasonable to assume that the temperature variation on a time scale r = 1 h cannot exceed 
1 K and 10 K on a time scale r = 12 h. Setting = (1 h) _1 = 1/3600 Hz in Eq. (p2|), we obtain then 7 = 3 and 
a = 1.41 K. 

Let us now estimate the precision required in the temperature monitoring in order that the amplitude of the 
temperature noise be lower than the "ordinary" noise ip(t). From Eq. (|20|), this requirement writes n 4 Av 2 6t(v) 2 < 
4>{vf. Taking into account that a{v) 2 = 8Q(v) 2 /Q 2 , = Cn 4 Av (Eq. (§)), and 7 = 3, one obtains from Eq. (|l]) 

the condition (for C = 2): 

fl <v^e^^=4.8xl0-f^f^) 3/ Y^V /2 K. (23) 

^ 300 KJ V 10" 5 Hz/ \ Ai< J v ' 

For v — (12 h) _1 , Av = 100 Hz, and O = 300 K, we get a < 1.7 x 10~ 2 K. This value is one hundred times smaller 
than the actual value of a obtained above from the expected daily temperature variation (a — 1.41 K). Since at this 
frequency ((12 h) -1 ), the amplitude of the temperature variation is about 10 K, this means that the temperature 
fluctuations have to be monitored with a precision of 10~ 2 x 10 K = 0.1 K for the extra-noise to be extracted from 
the "ordinary" noise ip(t). 



C. Periodic temperature fluctuations 



A periodic modulation due to the variation of the solar flux might be present in the power spectrum of the mirror 
temperature fluctuations. Maybe this will not be the case, given the amount of shielding the vacuum system and 
suspension will provide to the mirrors. We however consider temperature fluctuations correlated with the solar flux 
because they represent the most dangerous obstacle to the proposed method of detection. Indeed the solar flux 
spectrum contains, beside the solar day frequency, the sidereal day one, which may pollute the GW signal from 
galactic NS. Table || shows the Fourier spectrum of the solar flux at the latitude of Pisa (VIRGO site) computed by 
taking a period of 4 years. There are important lines at frequencies corresponding to multiples of the inverse of one 
solar day (f so i = 1.1574074 x 10 -5 Hz), as well as important lines at frequencies that are multiple of the inverse of 
one sidereal day (i^id = 1.1605763 x 10 -5 Hz). These latter lines can be explained by looking at the expression giving 
the solar flux as a function of time: the daily solar flux is modulated by the annual variation of the declination of 
the Sun. Consequently, there exists a line resulting from the beating of the solar frequency v so \ with the frequency 
corresponding to the inverse of the tropical year (vt y = 3.1688765 x 10~ 8 Hz). Therefore around the solar frequency, 
the two frequencies v so i ± v ty are present, and in particular the sidereal day frequency ^ S id = ^ so i + v t y 

The possible temperature fluctuation spectrum induced by the solar flux can be written 

5Q{ v f=Y J b 2 3 5{ V - V] ) , (24) 
i 

where the Uj are the solar flux frequencies displayed in Table O. 



6 



Let us compute the level at which the temperature of the mirrors has to be taken under control, in order to have 



the amplitude of this extra-noise lower than the "ordinary" noise ip(t). Taking into account Eq. (20), the relationship 
between the coefficients bj (Eq. (pi|)) and the periodic components aj of the noise a(t) reads 

(b 2 ) = (3 2 Q 2 &v 2 {a 3 ) 2 . (25) 

For a given multiple of the sidereal frequency, jv s id, we have thus to compare the contribution of the extra noise 
ajfi 2 {v) Ai/cos(27rjz4id£ + (/)) with the term ij}(t) (Eqs. (§) and (|l9l)). Taking into account that the thickness of 
the bin of a Fourier transform of a signal of length T is 2/T, we have a>j < ^20/ AuT; the corresponding periodic 
temperature variation yj (<50 2 ) must be less than 10~ 2 ( h z ) """^^yca/ 2 K (C = 2). 



D. Conclusions 



From the above analysis it appears that the non-periodic temperature fluctuations of the mirrors are not very 
dangerous, because an accuracy of 0.1 K in the mirror temperature measurement seems easily reachable. On the 
contrary, the periodic fluctuations must be kept under control within a few 10 -3 K for observation times longer than 
one year. This seems to be a challenging task. However, it must be noticed that the periodic temperature fluctuations 
can be measured on time intervals of the order of one year, and therefore are easier to control. Note also that measures 
of the noise at frequencies v ao \ and v ao \ — vty allows to deduce the noise at the sidereal frequency i^id = v so \ + vt y for 
the spectral components at the frequencies v so \ — z/ ty and ^ so i + v iy are almost identical (cf. Table ||). The required 
accuracy of these measurements is about 10~ 3 . 

To conclude, we recommend that all the detector parameters (temperature of the mirrors, power of the laser and 
so on) shall be monitored. In fact, as already said, we do not need any regulation of the temperature, but only to 
know, with the accuracy discussed above, how the temperature varies. 



V. NUMBER OF NEUTRON STARS AND THEIR GALACTIC DISTRIBUTION 



A. Number of rapidly rotating neutron stars in the Galaxy 



1. General considerations 



From the star formation rate and the outcome of supernova explosions, the total number of NS in the Galaxy is 
estimated to be about 10 9 |32|. The number of observed NS is much lower: ~ 700 NS are observed as radio pulsars 



150 as X-ray binaries [ p3|j34| (among which ~ 30 are X-ray pulsars) and a few as isolated NS, through their 
X-ray emission . 

For our purpose, the relevant number is given by the fraction of these ~ 10 9 NS which rotates sufficiently rapidly to 
emit gravitational waves in the frequency bandwidth of VIRGO-like detectors. The upper bound of this bandwidth 
(i/ max ~ a few kHz), is sufficiently high to encompass even the most rapidly rotating NS, at the centrifugal break-up 
limit, which — depending on the nuclear matter equation of state and on the NS mass — ranges from 1 kHz to 2 kHz 
J36|| . On the other hand, the lower bound of the interferometer bandwidth (v m in ~ 10 Hz) is a sensitive parameter for 
increasing the number of accessible NS. If the observed radio pulsars are representative of the population of rotating 
NS, lowering v m i n from 10 Hz to 5 Hz would increase the number of observable NS by a factor 2.3. 

In the following, we set the low frequency threshold of VIRGO to the value v mm = 5 Hz, which may be reachable 
in a second stage of the experiment. Using the fact that the highest gravitational frequency of a star which rotates at 
the frequency v is 2v, this means that the rotation period of a detectable NS must be lower than P max = 0.4 s. The 
NS that satisfy to this criterion can be divided into three classes: 

(CI) young pulsars which are still rapidly rotating (e.g. Crab or Vela pulsars); 

(C2) millisecond pulsars, which are thought to have been spun up by accretion when member of a close binary system 
(during this phase the system may appear as a low-mass X-ray binary); 

(C3) NS with P < 0.4 s but which do not exhibit the pulsar phenomenon. 
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In the examples given in Table |. the first three entries belong to class (CI), while the other two ones belong to 
class (C2). 

The number of millisecond pulsars in the Galaxy is estimated to be of the order Ns ~ 10 5 [3^] . The number of 
observed millisecond pulsars (P < 10 s) is about 50 and is continuously increasing. 

The number of young (non-recycled) rapidly rotating NS is more difficult to evaluate. An estimate can be obtained 
from the fact that the observed non-millisecond pulsars with P < 0.4 s represent 28 % of the number of cataloged 
pulsars and that the total number of active pulsars in the Galaxy is around 5 x 10 5 |3S[] . The number of rapidly 
rotating non-recycled pulsars obtained in this way is N± ~1.4x 10 5 . 

Adding N% and N2 gives a number of ~ 2 x 10 5 NS belonging to the populations (CI) and (C2) defined above. This 
number can be considered as a lower bound for the total number of NS with P < 0.4 s. The final figure depends on 
the amount of members of the population (C3). This latter number is (almost by definition !) unknown. We present 
below a scenario leading to a large and potentially detectable population (C3). 



2. The specific case of NS relics of the Galaxy formation 

It seems now well established that the birth of galaxies has been accompanied by a burst of massive-star formation. 
This took place at a cosmological redshift z ~ 2 [[59]], pQ ], i.e. when the Universe was ~ 20% of its present age. The 
remnant of these first generations of massive stars could contribute significantly to the population (C3), as we are 
going to see. 

The massive stars formed in the infancy of the Galaxy, some r ~ 10 Gyr ago, should have given birth to a large 
population of rapidly rotating (P _1 > 20 Hz) NS. Let us assume that these NS have a mean ellipticity of e ~ 10~ 6 , 
which amounts to only one thousandth of the maximum ellipticity permitted for the Crab pulsar (cf. Tabic |). 
Let us consider the fraction of these stars for which the spin-down is driven by gravitational radiation and not by 
electromagnetic processes. This may happen in the following cases: (i) the magnetic field is quite low (B < 10 10 G) 
either because the NS have been formed with such a field — given our poor knowledge of magnetohydrodynamical 
processes during stellar collapse and the proto-neutron star phase, this cannot be excluded — or because it has been 
destroyed during an accretion phase in a close binary system pH , (ii) the magnetic field configuration is such that 
the energy losses are small (for a review of the possible magnetic field structure of NS and their evolution see e.g. 
ref. [^U). Under these assumptions, the considered NS do not show up at present as radio pulsars, i.e. they belong to 
population (C3). Note that New et al. jl9| have also suggested that there could exist a large population of NS whose 
spin-down is driven by gravitational radiation. The differences with our hypothesis are that (i) they consider these 
NS to be presently rapidly rotating (millisecond periods) and (ii) they do not provide any specific scenario to create 
this population. 

An easy calculation show that with e ~ 10~ 6 and in 10 Gyr, the emission of gravitational radiation has slow 
down these NS to rotational period of P = 0.075 s, quite independently of their initial period. This corresponds to 
a frequency of 13 Hz, which fits in the low frequency part of interferometric detectors (contrary to the population 
suggested by New et al. |[9]]). The present gravitational wave amplitude emitted by such a NS at one distance unit 
is given by Eq. ( |A5| ) below and amounts to A ~ 7.5 x 10~ 28 kpc" 1 . Inserting this value into Eq. ( |l3| ) and taking 
the inverse-square average distance of this NS population to be D q = 5 kpc (cf. Sect. [II) leads to the following 
signal-to-noise ratio for the quadratic detection of these NS, the number of which is N: 

( N \ ( T \ /10 Hz \ 1/2 ( HT 21 Hz" 1 / 2 " * 



^ = 0.34 (_) (_)(_) ( ) , (26) 



where 10 -21 Hz -1 / 2 is the (present day) expected VIRGO sensitivity at the frequency of 10 Hz Q| and 10 10 is a 
possible value for the total number of relic NS. Note that since the considered NS population is supposed to radiate 
at low frequencies, one can take a narrow bandwidth Av — 10 Hz in order to increase the signal-to-noise ratio. Note 
also that improving the low-frequency detector sensitivity at from 10 -21 Hz -1 / 2 to say 10 -22 Hz -1 / 2 , would lead to 
a signal-to-noise ratio of ~ 30 if there are N ~ 10 10 relic NS with an mean ellipticity of 10 -6 . 



B. Galactic distribution 



The observed pulsars are concentrated toward the galactic plane, with a scale height above that plane of about 
0.5 kpc | p8[ . This latter value is almost an order of magnitude greater than the scale height of their progenitors 
(massive stars), reflecting the high "kick" velocities that pulsars acquire at their birth (see e.g. Q and references 
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therein). If the rapi dly ro tating NS considered in Sect. V A follow this distribution, their distribution function can be 
represented by Eq. (A21) and the corresponding time variation of the squared GW signal (h 2 {t)) is shown in Fig. [|. 
The Fourier spectrum of this signal is given in Table [II. Beside the constant part, it involves four frequencies, which 
are the first four harmonics of the sidereal day frequency v s id- The amplitude of the GW signal read in Fig. [j] is 



^(h 2 (t)} ~ 2 x 10" 




(27) 



This value can be compared with the GW amplitude of an individual NS which would have the same parameters as 
the mean ones used in the computation leading to Fig. 0, namely e = 1CP 8 , P = 5 ms and / = 10 38 kgm 2 [cf. Eq. (ffl)] 



h ~ 2 x 10 -27 



(28) 



Because of the important kick velocities mentioned above, it cannot be excluded that the actual distribution of NS 
constitutes a halo around our Galax y, instead of being concentrated in the galactic disk. A corresponding distribution 
function if then given by Eq. ( A22) and the resulting signal is shown in Fig. ||. 

From Figs, [I] and |, it appears that, from the detection point of view, the disk distribution is more favorable than 
the halo one. Indeed, the relative variation of (h 2 (t)} during one day is 100% in the case of the disk, but only ~ 15% 
for the halo. This is due to the fact that the halo distribution is more isotropic than the disk one. 



VI. CONCLUSION 



The detection of the gravitational waves emitted by rotating NS is a difficult task, but a positive result will 
give important informations on the evolution of massive stars. Two strategies are conceivable: the detection of the 
coherent signal from some isolated NS (linear detection) and the detection of the incoherent signal emitted by a large 
ensemble of NS (quadratic detection). The two strategies are complementary. The first one is well suited to detect the 
gravitational radiation emitted by a few very close NS, whereas the second one allows us to detect the gravitational 
radiation emitted by all the NS of the Galaxy. We have shown that if the distance of the closest NS is about 100 pc, 
the two methods give the same probability of detection if the number of radiating NS in the Galaxy is about 10 6 . 

The two strategies have their own drawbacks. Searching for individual radiating NS requires an extraordinary 
amount of computational time. It is possible that the overall sensitivity of this method will be limited by the power of 
future computers. The second method needs a high noise stability if it is to be used by only one detector. Moreover, 
this method has more chances to succeed if the number of emitters in the frequency bandwidth of the detector is 
large, which motivates the attempts to keep the low frequency threshold of interfcromctric detectors as low as possible. 
Besides, we have argued that a possible first generation of NS originating from an important stellar formation rate 
at the birth of the Galaxy, should have been spun down in 10 Gyr to periods of the order 0.1 s by the gravitational 
radiation reaction. This population could be detected by the quadratic method provided that the sensitivity of 
interferometric detectors is sufficiently good around 5 to 10 Hz. It should be noticed that since the wavelengths 
corresponding to these frequencies are between 3 x 10 4 km and 6 x 10 4 km, detectors at different locations onto the 
Earth can be employed for a search in cross-correlation. 

Finally, let us note that the quadratic technique presented in this article can also be applied to detect a single 
strong source of continuous wave gravitational radiation. 
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APPENDIX A: SQUARED SIGNAL FROM THE NS ENSEMBLE 

The detector output is 

s(t) = h(t) + n(t) , (Al) 

where n(t) is the detector's noise and hit) the detector's response to the gravitational radiation from the N galactic 
NS: 
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h(t) 



N 

£ 

1=1 



\F*(t)h\(t) + K(t)hUt)] 



(A2) 



In the above equation, h\ (t) and h % x (t) are the two polarization modes of the gravitational waves emitted by the i-th 
NS and F) (t) and F l x (t) are the (time- varying) beam-pattern factors taking into account the direction and polarization 



h l x can be 



of the radiation from the i-th NS with respect to the detector's arms (cf. Sect. 5.1 of ref. ||). h\ and 
expressed in terms of the the i-th NS's angular velocity f2j, ellipticity e^, distance r^, inclination of the rotation axis 
with respect to the line of sight tj and distortion angle Xii according to (cf. Eqs. (20), (21) and (25) of ref . [pi) 



A- 

h\{t) = — sinxz 

n 



Ai 



rl 



cosXi sin ii cos t, cos(flit + fa) 



sin Xi 



1 + COS Li 



cos2(f2ii + fa) 



rl 



•x\»;- -smxz 



cos Xi sin ij sin(r2jt + fa) 
sinxi cos ii sin 2(f2jt + fa) , 



(A3) 



(A4) 



A; 



c 4 P 2 



(A5) 



where Pi = 2tt /U,i is the rotation period of the star, Ii its moment of inertia, and fa some phase angle. 
Let us consider the mean value 



(h 2 (t)) 



1 



t+T 



h 2 {t') dt' 



of h 2 (t) on a time r such that 



£1 1 < t < 1 day 



(A6) 



(A7) 



where il is a typical NS rotation frequency: SI -1 < 1 s (for instance r = 10 s). From Eqs. (|A3])-( A4), and using (A7) 
(r 3> implies that integrals of products like cos(Oji) x cos(fiji) are negligibly small when i ^ j and t<1 day 
implies that F\. and F\ are approximatively constant on a time r) one obtains 



(h 2 (t))~h(t)+I 2 (t)+h(t) 



with 



' ^ 2 dt' 



h\ hi dt' 



JL (Pi \* rt+r 

<=1 T 174 



^ 2 dt' 



It can be seen easily that = -^(i) and, by virtue of (A.7), I/2I *C ii. Hence 



[/i l +(0] 2 dt 7 



(A8) 

(A9) 
(A10) 
(All) 

(A12) 



Using Eq. (AS) for h\ and performing the time integration leads to 

(h 2 (t)) — — [F l +{t)] 2 sin 2 Xi cos 2 x» sin 2 n cos 2 t 4 + sin 2 Xi(l + c °s 2 

i=i r * 



(A13) 
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For a given NS distribution, the expression ( A15 ) can be computed step by step, using the property that for two 
independent random variables u and v, and for large N, YliLi u i y i — ^Y2iLi v ii where u denotes the mean value of 
u: u ~ X)i=i u i I N . The variables A, \ an d <• being independent, Eq. ( |A13| ) results in 

N 

I ■ ■ I / 1 1 - 

(A14) 



A 2 



(h 2 (t)) = —— sin 2 xcos 2 \ sin 2 'cos 2 i + sin 4 \ (1 + cos 2 ii) 2 \ 



[F[(t)f 



Taking for \ and i distributions that correspond respectively to the probability law P{x) = 1/tt (uniform distribution 
in [0,7r]) and P(i) = 1/2 sint (uniform distribution on the celestial sphere) results in 



(A15) 



Fl_(t) depends on (i) the orientation of the detector's arms, which varies with time due to the Earth motion, (ii) the 
direction of the i-th. NS, which can be represented by its equatorial coordinates on the celestial sphere, namely the 
right ascension ccj and the declination Si and (iii) the polarization angle ipi of the gravitational wave with respect to 
the equatorial coordinates: 



(A16) 



The exact dependence is quite complicated and can be found in Sect. 5.1 of ref. g. Let us first take the mean value 
of Fl_(t) with respect to the polarization angle ip, which is uniformly distributed in [0, 2tt]. Let us denote the result 

by^F+: 



Equation ( Al5| ) then becomes 



+F+(ati,6i,t) = (F+(a l ,6 t ,tP,t)) i , 



/l2/ xx 43A 2 ^ 



240 



(A17) 



(A18) 



i=l 



Let f(r 7 a, 5) be the spatial distribution function of NS in the Galaxy, normalized so that N f(r,a,6) is the number 
density of NS: 



r—oo pS—Tr/2 pcx—2tt 



r=0 J 8=-iv /2 J a=0 



fir, a, 8) r 2 cos 8 dr dS da =1 



Equation ( A1S| ) becomes 
688 ir 4 G 2 



(h 2 (t)} 



15 c 8 



I 2 e 2 P- 4 N 



r—oo p8—tt/2 pa—27T 2 

/ / («, 5, t) / (r, a, 8) cos 8 drd8 da , 

r=0 J8=-ir/2 J a=0 



(A19) 



(A20) 



where use has been made of Eq. ( A5 ) to express A 2 . 

A Galactic disk distribution of NS can be modeled by the choice 



f(r,a,5) 



exp(-P/P ) 
AirR 2 , z 



if \z\ < z 
if \z\ > z Q 



(A21) 



where R — R(r,a,8) is the distance to the Galactic rotation axis and z — z(r,a,S) is the height above the galactic 
plane. Figure |l] shows the value of (h 2 (t)} computed according to this distribution with R = 3.8 kpc and z = 0.5 kpc. 
On the opposite, a NS distribution corresponding to a Galactic halo can be modeled by the choice 



f(r,a,S) 



G4 



(a/ao) 2 



(A22) 



637r 2 a 3 ^( ao / a ) _ l ' 

where a = a(r, a, 8) is the distance from the Galactic centre and ao is the radius of the halo. Figure || shows the value 
of (h 2 (t)) computed according to this distribution with ao = 10 kpc. 
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t [h] 

FI G. 1. Vari ation during one day of the total squared signal from a NS distribution concentrated in the Galactic disk, according to 
Eqs. (A.2C) and ( |A2l| ) with the parameters ~? = (10 -8 ) 2 , T 2 = (10 38 kgm 2 ) 2 , P~ 4 = (5 ms)" 4 , R = 3.8 kpc and z = 0.5 kpc. The 
localization and orientation of the detector are those of VIRGO, t = corresponds to a zero local sidereal time. 




t [h] 

FI G. 2. Var iatio n during one day of the total squared signal from a NS distribution corresponding to a Galactic halo, according to 
Eqs. (A20) and (A22) with the parameters ~? = (10 -8 ) 2 , T 2 = (10 38 kgm 2 ) 2 , P~ 4 = (5 ms) -4 , and a = 10 kpc. The localization and 
orientation of the detector are those of VIRGO, t = corresponds to a zero local sidereal time. 
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TABLE I. Gravitational radiation data for five selected pulsars. The GW amplitudes on Earth, ho, are computed according 
to Eq. (|l|) by assuming that / = 10 38 kgm 2 (a representative value for a 1.4 Mq neutron star). e„ 6 is the ellipticity in units of 
1CP 6 . The maximum ellipticity and maximum GW amplitudes are derived by attributing the totality of the observed pulsar 
spin-down rate to the emission of gravitation radiation. 



pulsar 


distance 




GW 


GW 


maximum 


maximum 




ellipticity 


name 




frequencies 


amplitude 


ellipticity 


GW amplitude 


to f 


;et ho = 10~ 26 




r [kpc] 


/ [Hz 


] 2/ [Hz] 


ho 


£max 








^detect 


Vela 


0.5 


11 


22 


1.1 x 10 _2Y e_ 6 


1.8 x 10 -3 


1.9 x 10 -24 


9.1 x 10" 


-6 


= 5 X 10 6max 


Crab 


2 


30 


60 


1.9 x 10- 27 e_ e 


7.5 x 10" 4 


1.4 x 10~ 24 


5.3 x 10" 


-6 


= 7 x 10"" 3 e max 


Geminga 


0.16 a 


4.2 


8.4 


4.7 x 10" 28 e_ 6 


2.3 x 10" 3 


1.1 x 10~ 24 


2.1 x 10" 


-5 


= 9 x 10" 3 e ma x 


B1957+20 


1.5 


621 


1242 


1.1 x 10 _24 e_ 6 


1.6 x 10 -9 


1.7 x lO -27 


9.1 x 10" 


-B 


^ Emax 


J0437-4715 b 


0.14 


174 


348 


9.1 x 10- 25 e_ 6 


2.9 x 10" 8 


2.6 x 10"" 26 


1.1 x 10" 


-8 


= 0.4e max 



a this is the recently determined distance from the measure of Geminga parallax 
b ref. || 



TABLE II. Fourier spectrum of the solar flux at the latitude of Pisa. 



frequency 


period 


identification 


spectral component 


v [Hz] ' 


u- 1 [h] 




(arbitary units) 


3.1688089E-08 


8766.000 


1 year 


0.6391195 


6.3376177E-08 


4383.000 


6 months 


1.4378662E-02 


1.1510698E-05 


24.13214 


^sol — 2z/ty 


4.0333651E-02 


1.1542386E-05 


24.06589 


Vsol — Vty 


0.4005030 


1.1574074E-05 


24.00000 


v so \ (solar day) 


1.836788 


1.1605762E-05 


23.93447 


^soi + fty (sidereal day) 


0.4006524 


1.1637450E-05 


23.86930 


faol + 2 fty 


3.9582606E-02 


2.3084773E-05 


12.03294 


2^sol — 2l<ty 


5.2575748E-02 


2.3116459E-05 


12.01645 


2f S ol — Vty 


1.2078404E-02 


2.3148148E-05 


12.00000 


2v BO i 


0.7087801 


2.3179837E-05 


11.98359 


2f so l + Vty 


1.4178075E-02 


2.3211524E-05 


11.96724 


2u so i + 2v ty 


5.2661020E-02 



TABLE III. Fourier spectrum of the total squared signal (h 2 (t)) corresponding to a disk distribution of NS and depicted in 
Fig, pj. 



frequency 


cosine coef. 


sine coef. 




(arbitary units) 


(arbitary units) 





7.895 





^sid 


0.923 


0.593 


2^ s id 


-0.738 


-0.065 




-0.079 


0.124 


4f si d 


0.005 


0.087 
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